NPS ARCHIVE 
1963 

CHIANG, S. 




i 




AN INVESTIGATION OP THE EFFECTS OF 
TRANSPORTATION LAG ON RELAY SERVO SYSTEMS 

SHU-MSN CHIANS 








„ VAl po X”oua. £ school 
CAUf0 ' N ' A 



AN INVESTIGATION OF THE EFFECTS OF TRANSPORTATION LAG 



CN RELAY SERVO SYSTEMS 



■>V 'k /V i: 



Shu-min Chiang 



AN INVESTIGATION OF THE EFFECTS OF TRANSPORTATION LAG 



ON RELAY SERVO SYSTEMS 



by 

Shu-min Chiang 

'/ 

Lieutenant, Chinese Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 



19 6 3 






\ \ 

n ieTC' 



LIBRARY 

US. NAVAL POSTGRADUATE SCHOOL 
montlrey, California 



AN INVESTIGATION OF THE EFFECTS OF TRANSPORTATION LAG 
ON RELAY SERVO SYSTEMS 
by 

Shu-min Chiang 



This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

/ 

ELECTRICAL ENGINEERING 
from the 



United States Naval Postgraduate School 



ABSTRACT 



Most physical systems are nonlinear to some extent. One variety 
of nonlinear system uses a relay or other on-off device as a control 
e lement . 

Since in the process industries or relay servo systems, one often 
encounters a type of time delay called transportation lag or dead time. 

It is highly desirable to investigate the effects of the transportation 
lag on the process or relay servo system. 

This paper is to demonstrate the effects of the transportation lag 
on the relay servo system and show how to eliminate these effects. 

Experimental results indicate that the transportation lag causes a 
continuous oscillation with a steady state over- and under-shoot, the 
response can never be settled, and also it leads to a longer rise time. 

In addition to those, the transportation lag results in a terminal limit 
cycle. Though the tachometer feedback has been used to compensate the 
system, this oscillation can never be eliminated completely when some delay 
exists. Therefore, a five percent overshoot and undershoot was considered 
acceptable . 

The author wishes to express his appreciation for the advice, assist- 
ance and encouragement given by Dr. Richard Carl Dorf, Professor of Electri- 
cal Engineering, United States Naval Postgraduate School, throughout this 
investigation. He also owes a debt of gratitude to Dr. John R. Ward, 
Professor of Electrical Engineering Department, U. S. Naval Postgraduate 
School, for his assistance during the initial stages of the work. 
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1. Introduction. 



Systems which utilize feedback for control purpose have become es- 
sential elements in modern technology. They range all the way from 
simple toys to the most complex automatic factories and production equip- 
ment. Feedback control is the guidance that makes modern automation pos- 
sible, and as such is in a large measure responsible for the ever-increas- 
ing productivity and rising standard of living of man. 

It is noted that feedback control systems are seldom truly linear, 
but rather they are subject to the introduction of a considerable number 
of nonlinearities. Some nonlinearities may be introduced intentionally by 
the system designer. An example of this might be the selection of a re- 
lay servomechanism with its attendent nonlinearity. On the other hand 
there may be no way of avoiding a nonlinearity in the system design. An 
example of this would be the mandatory inclusion of a gear train with its 
attendent backlash. A nonlinearity which frequently can neither be avoided 
nor neglected is tfansportat ion lag. 

In some control systems, especially in the field of process control, 

a non-eligible propagation time is encountered. That is, there is a definite 

lapse of time between the initiation of a change in the signal at one point 

in the system and the appearance at some other point of any effect of this 

change. But be sure that this kind of lag is not to be confused with the 

^ 2 

type of lag introduced by factors of the form 1/(TS + 1) or 1/ vT S +2 S TS 
+1) in the transfer function relating the two signals. Factors such as these 
in the transfer function of a component do reveal a tendency of the com- 
ponent to retard the signal, in a sense, by introducing phase lag at all 
frequencies, but they do not result in a finite dead time which must elapse 
before any effect of a change in input can be noted in the output of the 
component . 
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Dead time of appreciable magnitude can be observed, for example in 
long pipes carrying fluid under pressure, A change in pressure at one end 
of such a pipe is transmitted with finite velocity to the other end, pos- 
sibly but not necessarily being distorted in transmission. Another type 
of system in which dead time is often encountered is shown in Fig, 1-1, 

In such a system two different fluids are mixed at the mixing point, and 
the relative concentration is read at the metering point. If the meter- 
ing point is located at an appreciable distance from the mixing point, 
the dead time between the establishing of a change in relative concentra- 
tion and the detection of the change may be appreciable. 

Another example is illustrated below. In Fig, 1-2 the shaded area 
represents a slug of water at position 1, the outlet of the steam nozzle. 
Position 2 represents where the temperature sensing element will be placed. 
Let the temperature of the shaded area water slug at the steam nozzle be 
quickly changed from Q, to , Before the change can be measured, the 
slug of water must travel from the steam nozzle to position 2 where the 
sensing element would be located. 

The effect of the delay caused by the slug traveling down the pipe 
can be shown by locating sensing elements at position 1 and 2 as shown in 
Fig. 1-3, If the time constants of the sensing element at position 1 is 
taken as the input and the element at position 2 as the output, Fig. 1-4 
becomes the reaction curve for the process. 

This type of reaction curve is typical of a dead time lag (or trans- 
portation lag) and may be defined by the time lag T. The magnitude of 
the dead time lag will depend on the water velocity and the distance be- 
tween position 1 and 2, The lag time may be calculated as follows: 

T(min.) = d(ft) / V(ft/min) 
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Fig, 1-1 




Fig. 1-2 




input 

output 



t ime 






Fig. 1-4 
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where* T = dead time lag (min) 
d = length (ft) # 

V = velocity of flowing medium (ft/min) 

Dead time lag is very common in the process industries* Because 
of the complete lack of the response during the dead time* this type of 
lag is very detrimental to control. An examination of the equation de- 
fining the dead time lag shows that increasing the length or decreasing the 
velocity will increase the lag time* and we shall see later* make the 
problem of control more difficult. 

To determine the frequency response of the dead time lag* a sinusoidal 
wave having an amplitude of IQ V F will be generated around a bias tempera- 
ture 0^. At some low frequency the curve as recorded by the sensing ele- 
ments at position 1 and 2 would appear as shown in Fig. l-5(a). 

Fig. l~5(b) and 5(c) show similar results for two successively higher 
frequencies. As can be seen, as the frequency increases the period decreas- 
es. However* the lag time which is dependent only on the velocity and 
length remains constant. Thus the magnitude of the phase angle will in- 
crease proportionally to each decrease in period (increase in frequency). 

The phase angle may be calculated as follows: 

f> (degree) = terYodlScycle) X 360 (ddgree/cycle) 
or 0 (degree) = - T(min) X f(cpm) X 360 (degree/cycle) 

From this equation we see that the phase angle for the dead time lag 
is directly proportional to the frequency. The phase angle approaches 
zero degrees as the frequency approaches zero and negative infinity as 
the frequency approaches infinity. This will be discussed in jp 6. 
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Output Input Output Input Output Input 
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2. Relay servo system, 

2,1 General Description. 

When economy is of importance, many feedback control systems utilize 
a contactor for energizing the actuator. Such a system is called a contac- 
tor or, alternatively, an on-off control system, since the contactor is an 
element whose output can usually have only two values. In some cases 
these correspond to a fixed positive action and no actuator action. This 
usually is the case with thermostatic controls, for example. In other 
cases the two values correspond to maximum positive and negative actuator 
action. 

In Fig. 2-1 it is assumed that N is a relay or switch which has an 
output +1 or -I, depending on the sign of E = R - b. A positive output cor- 
responds to a positive sign of E and a negative output corresponds to a 
negative sign of E. When the switch is a perfect switch this can be de* 
noted symbolically by the equation 

V Sign E (2.1) 




Fig. 2-1 Block diagram cf the contactor servomechanism 
Since the contactor is frequently an electrically controlled or equivalent 
element, there may be a delay in the initiation of switching action as 
well as a finite time interval during which switching action takes place. 
The contactor action ordinarily applies full control force or torque to 
the system. Since the output of N is therefore not necessarily unity, the 



6 



actual value is included in the gain constant of the following element. 

The block G in Fig. 2-1 will be assumed to be a linear servomotor 
with viscous damping and no load torque except that due to inertia. The 
limiting case of zero viscous damping is of interest and will also be 
considered. Including the gain constant of the switch, the motor transfer 
function is 



G = 



1 Z 

ts( 3 t 



( 2 . 2 ) 



where T = KE is the stall torque of the motor and J is the total inertia 
o o 

of motor and load. In the absence of damping, 1 t'£ nj is set to zero in this 
equation. It should be emphasized that since N is nonlinear element, there 
is no particular reason to assume that the feedback element is linear. In 
particular, the feedback signal could be any algebraic function of position 
and its derivatives. However, for the method of analysis used here the 
feedback function can depend only on the position and velocity and must 
not on the derivative of the velocity (acceleration). The symbol H, there- 
fore, no longer represents a linear transfer function, but a possibility 
nonlinear feedback function involving the output position and velocity. 
Symbolically this will be written as 



b = f(C,C) 



(2.3) 



where C = dC/dt 
2.2 Relay servos 

The system behavior can be prescribed by combining the equations: 



o 

C 


+ 


/ 


* 

c 


c 


+- 


/ 

7~ 


c 



+ 



To 

J 

7Z 

J 



R -f( > o 

(2.4) 

R - f(c J c) <c 



Between switching times the response is described by a second order 
differential equation. For this reason such a system is called a second- 
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order contactor servo or feedback system. The two equations are inte- 
grable between switching times and the entire solution can be obtained 
by matching final and initial conditions of the two sets of solutions at 
the proper switching time. The solution is therefore easy to obtain in 
principle . The details, however, are so complex that another method of 
solution is desirable. 

The analysis of the system can be simplified by using the fact that 
it is possible to plot the response of the system as a trajectory in the 
phase plane, the trajectory being composed of a series of arcs or curves 
which have a characteristic shape dependent upon the properties of the 
system. The behavior of a physical system where only one coordinate is of 
interest is usually specified by determining this coordinate as a function 
of time, C = fj (t). This function depends upon the differential equation 
describing the system and upon a set of initial conditions. For some 
purpose, an equal suitable way of describing the behavior of the system 
is by means of a functional relationship between position and velocity: 




(2.5) 



When such a relation is known, the time dependence of C, if required, 
can be found directly by integrating the equation 



t/g 

cit 



c 



7 C (c.) 



Or 






rC 



c/^L 

f c v 



( 2 . 6 ) 



Equation 2.5 can be represented graphically with position and velocity 
taken as rectangular cartesian coordinates; position is usually plotted 
as abscissa and velocity as ordinate. The behavior of the system is des- 
cribed by a plot of all successive pairs of values of position and velo- 
city. The plane on which this plot is made is called the phase plane, and 
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the graph of Equation 2.5 is called the phase-plane trajectory. 

For the system to be considered in this paper, the differential 
equations of interest will not contain time explicitly. For the 
second-order system this implies that the acceleration can be written as 

c =F(c, c) 

The importance of this concept here results from the fact that for such 

systems, the behavior and therefore the trajectory are uniquely specified by 

# 

the initial values of C and C. This means that there is an only tra- 
jectory through a given (nonsingular) point in the phase plane. In this 
paper it is restricted to the systems described by second order differential 
equations and hence by phase plane trajectory in the two dimensional C-C 
phase-plane. 

The slope of the trajectory in the phase plane can be written as 

dc 

dc d t 

dc dc (2,7) 

dt 

According to this equation the slope is equal to the quotient of ac- 

«» 

celeration and velocity. In the upper half plane, C is positive and 
nonzero, and therefore the slope is nonifinite and either positive nega- 
tive, depending upon whether the acceleration is positive or negative. 
Similarly, in the lower half plane the slope is non-infinite and either 
positive or negative, depending on whether the acceleration is negative 

or positive. For a point on the horizontal axis the velocity is zero, 

* 

and therefore unless the accelerations also zero the slope dC/dC is in- 
finite. In general, then, trajectories cross the horizontal axis verti- 
cally. Direction of motion along the trajectory can always be determined 

e 

by the fact that for C positive C must become more positive with time and 
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for C negative C must become more negative with time. Motion along 
the trajectory is therefore to the right in the upper-half plane and to 
the left in the lower-half plane. The horizontal component of the 
trajectory velocity is equal to the system velocity C = dC/dt and 
the vertical component of the trajectory velocity is equal to the system 

c 

acceleration C = dC/dt. These considerations make it clear that a trajec 

tory can not end except on the horizontal axis. If there are points on 

* • •> 

the horizontal axis where both C and C are zero they are called 
singular points, and the value of C at such a point corresponds to an 
equilibrium condition of the system. Trajectories can emerge from or 
return to some of these singular points with a non- infinite slope, since 
the slope dC/dC is undefined at this point except through other considera 
tions, 

2.3 Phase plane trajectory. 

Some of the ideas in ^2,2 will now be illustrated with the system 
shown in Fig. 2.1 and defined by Equation 2.2 through 2.4, The Equation 
2,4, which are second-order differential equations, can be reduced to 
pairs if first-order equations by substitution of Equation 2.6 

k= C = r/C 

which gives 



Cft 



i/ f / i / / 7c 



( 2 . 8 ) 



an alternative expression for velocity, dV/dt = V(dV/dC), can be used 
with Equation 2,8 to give the following equation: 

V TFT = ~ i 



d 1 / , To 

^ = ~ zT ± ~r 



(2.9) 



This equation can now be integrated with the result 

c 7V 



v ± &Z LfO ( 2 . 10 ) 



10 



Evidently the constant of the integration is the value of C at 
which V = 0. This equation can be put in dimensionless form by the 
substitutions 

J 



1/ = V 

C„To 



C 



?jr„ 



c 



which gives 

z - z 0 



~[ v 2 ^ (7 f V) ] 



( 2 . 11 ) 



( 2 . 12 ) 



Typical trajectories for particular values of art shown in Fig, 2,2. 
Through use of the dimensionless variables v and c $ all systems of this 
type have identical trajectories. 

In the absence of damping, the trajectories are parabolas. If 
is set to zero, integration of Equation 2.8 gives the result 



c - Co =± 



X 



[/* 



2 



(2.13) 



In this case the substitution 2.11 can not be used. For the case of 
finite damping it is sometimes convenient to have an expression for the 
trajectory which shows the deviation from the parabolic shape more 
directly than Equation 2,10 or 2.12. This can be obtained by use of the 



power series expansion for log (1 + x) : 

To 



c _ c jL 



[It t 



TmTo 



_£ (JV_ 
4 * z,„ To 



A -7 



or 



c - 4 t ■§-, 1 / + 



<4 



o 2 

-(pjt 



(<?)*+■■■■ 7 



(2.14) 



(2.15) 



These equations are particularly useful when an approximate result is 
desired. For accurate computations Equation 2.10 or 2.12 is more use- 
ful, since the series given above converges rather slowly except for 
small values of v, and indeed diverges for jvj ^ 1. 

Reference to Fig. 2.1 and Equations 2.1 and 2,3 shows that the 
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Fig. 2.2, Phase plane trajectory for the second-order 
contactor servomechanism. 

condition governing switching from one trajectory family to the other is 
given by the equation 

/? - f (c, C) = o ( 2 . 16 ) 

For the present it will be convenient to restrict attention to step in- 
puts with initial velocity zero. In this case a simple translation of 
the origin can be used so that the initial conditions are c = -c , V = 0, 
and the final desired steady-state conditions are R = C = V » 0. Setting 
R = 0 in Equation 2.16 gives the equation 

( C * £ ) = O (2.17) 

The graphical plot of this equation separates the phase plane into two 
parts, in one of which the system follows trajectories corresponding to 
V - +1, and in the other of which the system follow trajectories corres- 
ponding to V = -1, this is illustrated in Fig. 2.3. 

For certain functions f(c,c) an anomalous situation arises. Consider, 
for example, the point P^ in Fig. 2,3, which will supposed to be reached 
by motion along the trajectory marked (1). When the switching line is 
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Fig. 2-3. Switching line for contactor servomechanism. 
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reached* the system will normally transfer to the trajectory of the 
other family which pass through P , this is the trajectory (2) in Fig. 

2.3. However* when the point is reached* for example along trajectory 
(3), a paradoxal situation arises. Although switching should take place 
at this time, there is no trajectory of the other family along which the 
system can move, since P^ is an end point rather than a starting point for 
trajectory (4). In terms of real physical systems the paradox can be 
resolved by recognizing that a small delaying switching is unavoidable. 
This permits the system to overshoot the switching line so that a small 
section of allowable trajectory is available in other family. When this 
trajectory reaches the switching line another slight delay allows over- 
shoot back into the left-half plane. The resultant trajectory consists 
of a series of short arcs in which the system oscillates about the switch- 
ing line as it approaches the equilibrium point at the origin. In terms 
of the ideal system, one can say that the trajectory follows the switching 
line to the origin accomplished by an infinite number of switching which 
provide, on the average, the proper force to give an acceleration consist- 
ent with the location in the phase plane and the slope of the switching 
line . 

Using the above considerations some interesting general conclusions 
can be drawn: 

(1) For no viscous damping ( — oa) and no rate feedback 

[•/Yc c) ^c] the trajectories are parabolas separated by the vertical line 

c = 0. The system is neutrally stable in the sense that for any amplitude 
of step input the oscillation will persist indefinitely, as shown in Fig, 
2.4a. 

(2) For viscous damping and no rate feedback j f the 
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Fi &° 2 ' 4 * Typical phase plane trajectories 
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trajectory are the "asymmetric parabolas" of Fig. 2.2, separated by the 
verticle line c = 0. The system exhibits a damping or decrement of the 
oscillations produced by a step input. The per-unit decrement is a func- 
tion of amplitude and approaches zero as the amplitude approaches zero. 

In other words, for large input the initial per-unit decrement is high, but 
decreases to zero as the amplitude of the oscillation approaches zero. 

This can be observed qualitatively in Fig. 2.4b. 

(3) With linear rate feedback, the feedback function is ffc* 

If there is no viscous damping ( , then the trajectories are 

parabolas separated by the straight line 

C = - -fi- C (2.18) 

Again the system shows a per-unit decrement or damping which varies with 
amplitude. However, in this case the damping is small for large amplitude 
and as the amplitude decreases the damping increases until the final tra- 
jectory into the origin lies along the line 




with no further overshoots. An illustration of this is shown in 2.4c. 

The above conclusions shows that no simple linear feedback function 
will give optimum response, regardless of the system's viscous damping. 

If rate feedback is adjusted to give no overshoot for the largest expected 
input then the system moves to the origin slowly for small inputs. If the 
rate feedback is adjusted to give faster response for small input, then 
overshoot and oscillations will follow application of a large step input. 

It is natural to ask whether any feedback function can be found to 
give optimum response for all inputs. This turns out to be a difficult 
question. Since the system is nonlinear and does not have an easily deter- 
mined response for arbitary inputs. However, a definite answer can be found 
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for the case of step inputs by following chain of reasoning. For most 
rapid response without overshoot it is necessary to apply full positive 
acceleration torque until that instant at which full negative torque will 
bring the system to rest at the origin in the phase plane. This is il- 
lustrated in Fig. 2,5 for a particular initial starting point. The tra- 
jectory shown passing through the origin can be designated as the homing 
trajectory. In order to achieve switching on this trajectory for all in- 
puts, it is necessary that feedback function set equal to zero, 

j (<z, a) — O (2.19) 

being the equation of this homing trajectory. In this case, when any 
trajectory crosses this line switching will occur, so that the applica- 
tion of full reverse torque will just bring the system to reset at the 
origin. 




Fig. 2.5 Switching point for optimum step response 
Since the trajectory of Equation 2.8 is roughly parabolic, a suit- 
able form of the feedback function is 

'f C c ' v) ^ k ( C t \ )= C ~+ V ji/| 
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Strictly speaking, this equation should be modified to match the appropri- 
ate "asmmetric parabola”, although for practical purposes the refinement is 
not necessary. In other words, if the position feedback is linear and the 
rate feedback roughly quadratic (except for sign) then it is possible to 
obtain optimum response for this particular type of input. It will be 
seen that switching delays give rise to oscillatory motion of small ampli- 
tude unless a linear velocity-feedback term is introduced. For this reason 
it is usually desirable to introduce a small amount of such feedback (in 
addition to the quadratic feedback) in order to eliminate or minimize this 
oscillation. This does not significantly affect the results for large 
amplitudes . 

It is quite important to remember that this optimization does not 
necessarily give optimum response for arbitary inputs. Further considera- 
tion is beyond the intended scope of this paper. 

2.4 On-off control. 

There probably are more on-off discontinuous controllers in the world 
than the number of all other servomechanisms and regulators put together. 
Hot-water heaters, ovens, refrigerators, furnaces, electric irons, battery 
chargers, liquid-level controls, and compressed-air pressure regulators 
are all commonly of this type. 

A relay system is one which contains a decision function. It is 
commonly in the error amplifier and in its most elementary form is an on- 
off switch. If the positive error exceeds a certain amount, the switch is 
turned from on to off, and if the negative error exceeds a certain amount, 
the swtich is turned from off to on. When the switch is on, a predetermined 
constant corrective force is applied into the rest of the system to bring 
the output back into correspondence with the desired value. When the switch 
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is off, this force is removed. The differential equation for each is the 
open-loop equation, with one term, the corrective force, having two differ- 
ent values for the two periods. Each switch operation starts a new period 
for which the initial values are the final values of the previous period. 

One characteristic of a relay is that it delivers an output having one 
of only two values, independent of the continuous changes in the input, 
except they specify the switching instants. It also has hysteresis, the 
phenomenon that the input values for switching on and for switching off 
are not exactly equal. The actual time of motion of the armature will be 
considered negligible, but the system may contain the dead time of flow of 
material through a long pipe, equivalent to the time of transmission of an 
electric signal on a long distoritonless line. 

The term relay is used to apply to any device having at least the first 
of these characteristics. In this paper, we are interested primarily in 
discontinuous control of continuous processes, so that the relay is only 
one part of the system. 

The on-off control, or two position control, is undoubtedly the most 
widely used type of control for both industrial and domestic service. It 
is the kind of control generally employed on home heating systems and 
domestic water heaters. 

Two position control is a position type of controller action in which 
the manupulated variable is quickly changed to either a maximum minimum 
value depending upon whether the controlled variable is greater or less 
than the set point. The minimum value of the manupulated variable is usually 
zero (off). This mode of control is illustrated by the electric level- 
control in Fig. 2,6. A float in the vessel operates an electric switch 
which controls power to a solenoid valve. When the liquid level rises, 
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Fig. 2-6. Two position control 
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the switch contacts are closed, the solenoid valve closes, and the flow is 
cut off. When the liquid level falls, the switch contacts are opened, the 
solenoid valve opens, and the inflow resumes. If the float lever has no 
bearing friction and the electrical contact draw no arc, the action is sharp 
or "knife 1 edge” as shown in Fig. 2-7. 

The equation for the two-position control are 

m = when e P'" 0 

m = M when e 0 
o 

where m = manupulated variable 

= maximum value of manupulated variable (on) 

= minimum value of manupulated variable (off) 
e = deviation 

Thus two-position control must be described by two equations, each apply- 
ing in a certain region of deviation. 

A differential gap in two-position control causes the manupulated vari- 
able to maintain its previous value until the controlled variable has moved 
slightly beyond the set point. In actual cperation it is the same as hystere- 
sis, as may be seen in Fig. 2-8. 

A differential gap is caused in the two-position controller of Fig. 2.6. 
If small static friction exists at the bearing on the float arm, the liquid 
level must rise slightly above the desire value to create sufficient buoyant 
force to overcome friction when the level is rising. Also, the liquid level 
must fall slightly below the desired value when the level is falling so that 
the weight force may overcome the friction. This kind of differential gap 
may be caused by unintentional friction and lost motion. 

A differential gap may be intentional, as when a magnet is installed on 
the float arm in Fig. 2.6, causing a hysteresis in float-arm action. Similar 
arrangements are common in domestic thermostats and are employed for the 
purpose of preventing rapid operation of switches and solenoid valves for 
reducing of electrical contacts. 



21 



3 . Description of the system . 



The nonlinear feedback control system to be studied is a relay servo 
system as shown in Fig* 3-1. 

Block one contains the transportation lag. The nature of this block 
and its analytic representation will be discussed in section 6. 

Block two is the ideal relay, no dead zone and no hysteresis* 

Block three is the motor load unit or plant. It is assumed to be a 
second order system* In order to have the transfer function always low 
pass in nature, we choose the value of G(s) as 



Block four is a tachometer. In this system velocity feedback compen- 
sation is used. 

This system, in addition to the velocity feedback, will also use unity 
position feedback. 




/ 




(3.1) 
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Fig. 3-1 Generalized block diagram 
of basic system 



4. 



Nonlinear analysis . 



4.1 General 

In many cases the system designer is forced to work with components 
which are not entirely linear. In order to linearize the response of a 
nonlinear actuator we use the feedback. The basic difficulty with the non- 
linear equation is that it is not possible to exhibit a general solution. 
Associated with this is the fact that there is no complementary function 
or transient solution made up of terms with characteristic shapes, there 
are, therefore, no system roots such as were defined for linear systems 
with constant coefficients. It is not the scope of this paper to include 
a comprehensive treatment of nonlinear systems. 

Many systems are approximately linear when the amplitude of the 
system variables are not too large or too small. For large enough ampli- 
tude all physical systems must become nonlinear. A system which for small 
amplitude is a linear unstable system must become nonlinear when the oscil- 
latory amplitude exceeds some value. Such a system will sometimes reach a 
stable oscillatory motion which is termed a limit cycle. Although such 
motion is usually periodic, it need not be sinusoidal. 

There are other large-amplitude nonlinearities which are in a sense 
destablizing. In these systems, small amplitude response may be stable 
and satisfactory, but large disturbances may send the system into uncontrol- 
led large-amplitude oscillations. 

There are also nonlinearities which become important for small amplitude 
motions. Mechanical friction and backlash are good examples of this type. 

In the case of backlash it is not unusual to find a limit cycle of small 
amplitude which leaves the system capable of following commands or inputs 
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with good accuracy and in an essentially linear mannes except for the small 
limit cycle motion which is superimposed on the desired response. 

4.2 The describing function. 

Many functions involved in automatic control are not desirable by 
elementary techniques, and they must be described by graphical or other 
methods. 

The output of a nonlinear element will not necessarily be sinusoidal 
even though the input is a pure sinusoid, but it will generally contain 
fundamental and harmonic waves. The frequency of the fundamental of the 
output will often be the same as the frequency of the input. In this case 
a describing function can be used and the element may be treated as a lin- 
ear element. 

Since most elements of a feedback system are of a low-pass nature, the 
higher harmonics may be greatly attenuated in the rest of the loop, so that 
the input to the nonlinear element may be quite accurately sinusoidal. If 
this is the case, it is necessary to compute the transmission around the 
loop for the harmonics. In particular, for the nonlinear element we need 
only a function to describe the fundamental component in the output as a 
function of the amplitude and frequency of the sinusoidal input. This furc- 
tion has been named a describing function. Since the describing functions 
are approximations to the output-input ratio of an element having slight 
nonlinearily, it is anolagous to the transfer function of a linear system, 
and for certain purposes will be used in a similar manner. There is, of course, 
one important difference. Since a transfer function describes a linear sys- 
tem, it may be frequency-dependent but it is not amplitude-dependent. On 
the other hand, the describing function is usually amplitude-dependent and 
may also be frequency-dependent. The describing function is expressed in 
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terms of the amplitude and phase of the fundamental of the output signal 
related to the amplitude and phase of the input signal. The describing 
function can be obtained analytically or experimentally. 

Describing functions are generally useful in automatic control anal- 
ysis if the neglect of the harmonics in the output signal does not cause 
difficulties. They may be used with adequate results when the element in 
question is followed by one or more elements having appreciable attenua- 
tion of the harmonics in question. 

4.3 Computation of the describing function. 

Consider a nonlinear element whose input is sinusoidal, 



x = X sin Oft 



(4.1) 



and whose steady state output can be written as a Fourier series: 




b^ sino/t + a^coslc^Jt + b^ sin 2aJt 



+ 



= a I 2 — ^ (a cosk^ut + b sink^t ) + ....(4.2) 

O , t- K K 



Evidently 




(4.3) 



b 




(4.4) 



The describing function is defined as the complex number 



N = Nj + jN 2 = N t € 



jf 



(4.5) 



where 




N 2 = a x /x 



or, alternatively, 




(4.6) 
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This definition immediately puts certain restrictions on the permissible 
nonlinearities. The output function must be a periodic function with 
period equal to that of the input. According to this definition, devices 
which produce subharmonic response do not have a describing function. 

The form of Equation 4.1 implies that the input is sinusoidal and 
does not contain d-c bias. If a nonlinearity in a feedback system is not 
symmetric, then in general there will be such a bias. 

Because we are only interested in the nonlinear device of ideal relay 
only ideal relay will be discussed. The characteristic of the ideal relay 
is as shown in Fig. 4-1. 



+ V 
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Fig. 4-1. Characteristic of ideal relay. 

From Fig. 4-1 



e = + Y o ujt ^ it 

o 

e = - Y rt <^t< 2 re 

o 

The describing function of the ideal relay can be obtained very simply 
by carrying the integration indicated in Equations 4.3 and 4.4 and 4.6. And 
in calculating the describing function, we are only interested in the funda- 
mental terms. Thus we have the describing function of the ideal relay as 



QJ/oj)= 



4 y 

7TX 



(4.7) 
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From the above results, we know that the describing function of the 
ideal relay is amplitude-dependent and independent of frequency. This 
function is plotted against X in Fig. 4.2. 




Fig. 4.2. Describing function of ideal relay. 

4.4 Analysis of the limit cycle. 

4.4.1 Frequency response method 

Limit cycle angular frequency , as predicted by the frequency re- 
sponse method, is such that the loop phase shift is 360 degrees and this 
occurs when the loop frequency response function has no imaginary part. 

The flow diagram for sinusoidal quantities is shown in Fig. 4.3; the 
relay is represented by a describing function in cascade with a time de- 
lay transfer function and the relay output takes on the values + or - unity. 
The condition for 




Fig. 4.3. 

a limit cycle is therefore met when 



Q 



-Jos*, td 



is real 



</*£ (ju!± -+ /J 

by rationalizing this expression the condition becomes: 
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is real 



Expanding the exponential term gives the limit cycle condition as 



(Cos - jSin C<J Z t d )(6jJ + /“U is real 



Equating the imaginary terms to zero yields the condition 



OJl fee* COi = / 



P 



2 7T 



(4.8) 



2n 



For a given value of t , (delay time) the dimensionless period P = —— — 

d 

is multivalued. These values correspond to the possible modes of oscilla- 
tion and only the longest period is stable. 

According to the frequency response method the servo loop gain at limit 
cycle frequency is unity since the oscillation is exactly regenerated in the 
loop. To determine the limit cycle amplitude at the input to the on-off 
element it is therefore necessary to solve the equation obtained by equating 
the loop frequency response function to unity at the angular frequency, that 



IS 



2 « j 

n X u l JoJ L - OJj? J 1 



(4.9) 



where is the amplitude of the fundamental component of the error signal 
when limit cycling. Thus at limit cycle frequency, the imaginary part of 
the loop frequency response function is zero and by rationalizing the left- 
hand side of equation 4.9 and equating the imaginary part to zero we obtain 
4 j C-Q$ td -t oJj. 

l 



4 

7TX, 



or 



/ ■+ 






X = Y J=r -fL p J && Cos tci L 7 

A/. ' J r+ co-j / < 4 - 10 > 

where x^ is the dimensionless measure of misalignment limit cycle amplitude 

and P is the time limit cycle period. 

The values of P and x from equations 4.8 and 4.10 are plotted as shown 

L 
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in Fig. 4,4. From Fig. 4.4 we can predict limit cycle amplitude and period 
in an on-off system. 

Equation 4.9 can be rewritten as 



G ^ ( 4 - 11 ) 

Equation 4.11 may be analyzed for stability by Nyquies method or by root 
locus method if the describing function is amplitude-dependent only. 

When is amplitude- but not frequency-dependent, it is possible to 
reduce a very simple criterion for predicting stability or instability. 

The basis for the criterion is that the leftside of equation 4.11 is only 
frequency-dependent and the right side is only amplitude-dependent. Margin- 
al stability as predicted from the Nyquist criterion and the describing fun- 
ction concept requires a solution of equation 4,11, or in other words, an 
intersection of these two curves. At the point of intersection, the frequency 
parameter on the one curve and the amplitude parameter on the other give the 
frequency and amplitude of the corresponding oscillation. If there is no 
intersection, then the system is presumably stable or unstable for all condi- 
tions of operation. Since it is difficult to find a rigorous definition of 
stability for nonlinear systems, the above statement must be regarded as pro- 
visional, pending more detailed investigation or more careful definition of 
stability for urusual systems. 

From equation 4.11, we have ■ '** — 

j vj ( -r t) 



/ 7T X . _ __ . 

Qo 4-Y /' 5 70 $ X 

if we put t^ = .2 seconds and by assigning some values of zxT , we can have 

the following data; and from these data the Nyquist plot can be plotted as 

shown in Fig. 4.5. 
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Fig, 4-4 Curves for predicting limit cycle amplitude and 

period in r . n on-off system 
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From the Nyquist plot, a limit cycle is observed; it is at a point where 
,52 and { 1*71. (The intersection of these two curves). The 

intersection may indicate a stability boundary in the sense that excitation 
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at level below that corresponding to the intersection will be unstable for 
the system, and should reach a limit cycle in which the input to the non- 
linear element is roughly sinusoidal with angular frequency OJ ~ * 52 and 
X = 1.71. The system should be unstable for low level of excitation and 
stable for level above that corresponding to the intersection. 

4.4.2 The root locus method 

Because the describing function of an ideal relay is independent 
of frequency, so that the root locus method can be readily applied. The 
nonlinearity introduces no phase shift but only a gain dependent on ampli- 
tude. It will be assumed here that the behavior of the nonlinear system 
can be approximated through use of a root-locus diagram for the system with 
the nonlinearity eliminated. By block diagram manupulation as follows; 




( a ) 




Fig. 4.6. 

and algebraic manupulation provides the stability relationship, the 
characteristic equation 

Gj G 2 G d (H +1) = - 1 (4-12) 
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G(S) = - 1 



(4.13) 



G(s) = K ' V 4 e ' St 
S(S + 1 ) 



k' = 



4Y 

nX 



For t = .5 second, 
d 

put K t 



K fc S = H 



=1.5 then we have 



_ K (S + 0.67)e 
G(S) ~ S (S + 1) 



- . 5S 



Equation 4.13 is complex and may be split into two equations by equating 
the magnitude and phase angle on both sides respectively. Thus we have the 
magnitude equation 

| G(S)| = 1 (4.14) 

and the phase angle equation 

/G(s) = ± 180° (4.15) 

where 



where 



/G(s) = /g'(s) ~ Ctrl 

/g'(s) -Ml = + 180° (4.16) 

g ' (s, = ■ * 

V ' S (s + 1) 

a method is introduced to obtain the required locus of equation (4.16). 
First it is to construct a family of root loci (or phase-angle loci) of 
the equation 



SS ' (s) = 0 X 

where 0^ is a constant phase angle. Thus a family of loci can be plotted 
by various values of 0^, Fig. 4.7 illustrates a family of root loci for the 
transfer function G t ( c; ). Since the locus with phase angle 0, is symmetrical 



34 



b 




35 



Fig. 4-7 Construction of the Required locus with 

Transportation lag 



to that with phase angle (-0) with respect to the real axis, only the loci 
on the upper half plane are shown. 

Next, it is to construct another family of phase-angle loci (the name 
dead time lag loci is suggested) on the s-plane from the relation 

-CO T = 0 2 

the dead time lag loci family is also plotted in the same sheet as shown in 
Fig. 4.7. Then by superimposing these two families of loci on the same 
s-plane, the point of intersection of the corresponding curves with the sum 
of phase angles 0^ and 0^ equal to 180 degrees are the points of the locus 
required. The curve drawn through these points is the required locus. It 
is the required locus because each point on it satisfies equation 4.16. 

The determination of the gain can be carried out then by the criterion in 
the root- locus method. The magnitude equation is 



K " |CH<r + j of ) | e" Ts (4.17) 

which will give the gain K of any point on the required locus. 

The required roots are the points where the value of gain K is equal 
to one. The gain limit, which is the maximum gain below which the system 
can have damped oscillation, can be found readily. This is the value of 
gain at the point of intersection of the required locus and the imaginary 
axis. This value is found to be 1.21. 

This given function transfer has an infinite number of roots, which 
can be obtained by constructing additional branches of required locus. 

In general, the number of influential roots depends on the relative 
magnitudes of the dead time T and the other time constant in G'(s). 

The limit cycle occurs at CL/ - 1.67 rad/sec. and k = 1.21 where 



k = k f k 



K 1 = 



4Y 
n X 



2 _ 

JT X 
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k = 1.5 



K = 1.21 = 
X = 



n X 
2 x 1.5 



« x 1.21 



K t n X 
= .788 



x 1.5 
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5. System with Dead Time , 

5*1, Transfer function associated with pure dead time 

If the input to component with a pure dead time characteristic is f(t) 
u(t), then the output of the component is f(t - T), where T 0 is the 
dead time/ The transfer function associated with pure dead time is, there- 
fore, 



G(s) - 



-TS , N 
e F(s) 

F(s) 

-TS 



For sinusoidal input at frequency clt 
in equation 5.2 becomes 

i 

r» 



(5.1) 

(5.2) 

the transfer function given 



G(j ur) = (5.3) 

The gain is thus seen to be independent of the frequency and equal to unity, 
and the phase shift is seen to be directly proportional to frequency. 

5.2 Frequency response of a system with dead time. 

It is of interest to compare the frequency characteristics of a unit 
with pure dead time T and a unit with a first-order lag transfer function 
with time constant equal to T. That is, 

-TC 

G x (s) = e (5.4) 



and 

c 2< s > = tsVt <5 ' 5) 

The gain associated with G^CjcoO and GjCj^cO are 

| G l ( j^)| = 1 (5.6) 

r ■ 5 n "j 

j G 2 (j os )| = V 1 + cal (5.7) 



and their respective phase angles are 



0 l = -a/ T 



(5.8) 
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Fig. 5-1 Bode diagram for a pare dead-cine unit and a ?i. rec- 
order lag unit with time cor.r.t, i T :d: Vtto 

equal to the dead time of tr 
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and 



0 ? = - arc tan OJT (5.9) 

The Bode diagram for the two units are plotted as shown in Fig. 5.1, where 
it can be seen that for a?T 1, the units have nearly the same effect on 
a given input signal. Generally, however, since for higher frequencies the 
dead-time unit exhibits more phase lag and more gain than the first-order 
lag unit, the introduction of the dead-time unit in cascade within a given 
closed-loop system has more adverse an effect on the stability of the system 
than does the introduction of the first-order unit. 

Dead time has the characteristic of a pure phase shift since the magni- 
tude is always one. Seen easily from the graphs drawn in Fig. 5.2 that no 
other dynamic element builds up phase lag so much as dead time. In combina- 
tion with other dynamic elements the phase lag of dead time is simply added 
to obtain the total phase lag. 





Fig. 5-2. Dead-time element. 

A useful approximation for dead time T is the following function from 
the Pade list, 



, f os _ 2 - TS 

G (S) - 2 + TS 



(5.10) 



40 



The magnitude of this function is everywhere 1.0 and the phase lag cor- 
responds very closely to dead time up to frequencies given by a/I = 0.5. 

In addition, this function may be easier to incorporate into complex system 
analysis . 

The destabilizing effect of the transportation lag, as dead time is 
also known, and clearly illustrated in a Nyquist diagram. For example, if 
a unit with a pure dead-time characteristic given by 

Gj (S) = e' TS (5.11) 

is connected in cascade with a component with transfer function 

y s > - s(s 2 +n (5 - i2 > 

as in Fig. 5.3, the Nyquist diagram for the system for various values of 
the dead-time appears as in Fig. 5.4. As can be seen in this graph, for 
small values of transportation lag the introduction of the dead-time unit 
merely results in a decrease in relative stability, but for large values of 
transportation lag the system can be unstable. 

The effect upon stability of the introduction of a transportation lag 
in cascade within a closed-loop system of the type illustrated in Fig. 5-5 
(a) can be determined from the Nyquist curve for the original system in the 
following manner. The characteristic equation of the original system is 

1 + G a (8)G b (s) = ° (5.13) 

and that of the system which is obtained by introducing a transportation 

lag T within either the forward path or the feedback path is 

1 + G a ( s ) G b (s)e " TS = 0 (5.14) 

Equation (5.14) can be rewritten as 

G a (s)G b (s) = -e +Ts (5.15) 

or, in terms of frequency response, as 

G a(j^) G b (j ^) =*e j " /T (5.16) 
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R + 



®r-~ 







Fig* 5~3o Block diagram of a simple feedback 
system with transportation lag. 




Fig, 5-4, Nyquist diagram for the system 
illustrated in Fig* 5-3, 
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( a ) (b) 



Fig. 5-5 Block diagram of a single loop feedback system 




Fig. 5-6 Stability locus in the G G -plane for a 

Q. D 

system of the type shown in Fig. 5-5(b). 




Fig. 5-7 Determination of phase margin in a 
system with transportation lag. 
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The locus of G (jcc/')G, ( j ckT ) is the Nyquist curve for the original system; 

a b 

equation (5-16) states that the effect of introducing transportation lag T 

is to shift the critical stability point in the G G, -plane from -1 + jO to 

a d 

-coso/T - jsimxrT. That is, the critical point becomes a critical locus, 

in the form of the unit cycle in the G^G^ -plane, as shown in Fig. 5-6. 

The characteristic equation (5-16) is satisfied only if the locus of 

G ( jc</)G, ( j cu ? ) intersects the stability locus in such a manner that the 
a b 

frequency on the Nyquist curve and the frequency on the stability locus are 
the same at the point of intersection. If the two frequencies are not equal 
at the point of intersection, a measure of relative stability is afforded by 
the stability locus. Thus, if the Nyquist diagram for a given system has the 
appearance shown in Fig. 5-7, the phase margin can be determined as follows: 

1. Determine the cross-over frequency , at whichthe Nyquist 

curve crosses the unit circle. 

2. Locate the point on the stability locus for ccr = clS^ . 

3. Draw a radial line segment from the origin through each of 
these points. 

4. Measure the angle between these line segments. This angle is 
the phase margin, with the same positive sense as in the case of a system 
without transportation lag. 
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6. Transportation lag . 



6.1. General 

The phenomenon of dead time, sometimes called transportation lag or 
time delay or ideal time lag, can best be categorized as a lack of change 
in output for a specified interval of time following a change in input, 
this is illustrated in Fig. 6.1. 




Fig. 6.1 

Transportation lag may be encountered in process involving pressure, 
temperature, and also may be found in process involving chemical reactions 
when a finite time must elapse before a reaction begins to occur. 

Occasionally in analyzing engineering systems, it becomes necessary to 
simulate time delay (or transportation lag). For example in studying trans- 
mission processes through hydraulic or pneumatic lines or in studying the 
motion of radio active particles in neuclear reactors it is necessary to 
produce an output voltage which is identical with an input voltage but dis- 
placed in time by a specified amount. 

Transportation lag may be expressed by the following mathematical re- 
lationship: Given f(t) for t ^rO, the function f(t - T) differs from the 

function f(t) by the dead time T, where f(t - T) is defined to be zero for 
t < T. 
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A direct and accurate method for generating such time delay is to 
employ a magnetic tape recorder in which the space between the recording 
and reading heads is adjustable, by suitably controlling the relative loca- 
tion of the two heads and the speed of the tape a wide range of time delay 
can be simulated. 

Approximate methods which do not require special auxiliary equipment, 
once an analog computer is available, are based upon the Laplace-transforma- 
tion relationship 

/[ / (t - Tj] = Ffs j e 5T = U[ f«>] e' r (6.1) 



- sT 

since e cannot be expressed as a rational function, the transfer func- 
tion for dead time can be simulated only approximately, thus the remaining 

c 

methods depend on an approximate series or technique of network synthesis,"' 
The Taylor series expansion of this function is 






/ - TS -t 



(TSX 
2 I 



(TS) d , (Tsf 
3 1 ' 4 / 



f 



(TJ )' 1 

a / 



( 6 . 2 ) 



The generation of the terms of this series would require successive differ- 
entiation with attendent noise problems. The rate of convergence of equa- 
tion 6,2 is further relatively slow, particularly for large values of T, 
or high frequencies. For this reason the Taylor series is not well suited 
to the generation of transportation lag for problems involving high fre- 
quencies or long time constants. 

6.2 Pade approximation 

A more practical approach to the generation of the transportation lag 

9 

function is to employ the so called Pade 1 approximation. From equation 
6.1 the transfer function of transportation lag is 

(C =f ST 'A ^2. (6.3) 
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substituting s = jar to find the frequency response gives 



e JUJT = Cos«jT - jS.nco'T 

E c 



thus 



ft 



= Jc&S 2 as T -f S>'n 3 co / 



■= fan '(- 



5 is? oo_ 7~ 

Cos as T 



■) 



■ / 

— co ~ 



(6.4) 

(6.5) 



from the above equations the frequency transfer function of the trans- 
portation lag is characterized by an amplitude ratio which is unity for all 
frequencies and a phase angle which varies linearly with frequency over the 
desired range of interest. From equations 6.4 and 6.5 it is seen that the 
quotient of any pair of complex-conjugate polynominals has an amplitude of 
unity and also possess considerable phase shift. 

In the Pade' approximation the exponential is approximated as the ratio 
of two polynominals in s, according to 



P„jS) 



( 6 . 6 ) 



where subscripts m and n refer to highest power of s in the numerator and 
denominator, respectively. Expressed in its basic form, the Pade* approxi- 
mation is 



<S? X Li'm (2Ll 

(o+v) — oo Gru tV (x.) 



where 




l/X V(y-QX*_ 

(O-H/) 1 . (<-J + v)(o-h/->)2 / 






V(v-_ 0(w-2) - - • . • / X ^ 

(<-/+ i/) (o+ is-/) (CJ + I ) |/ I 



(6.7) 



( 6 . 8 ) 
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o(u-/) 



/ (if to)! ^ (V+U) (ift- 

r ,\V u (Q-oco- 2 ) *?- / X * 

*+ ( v (y+u )(i/+u- 4) ■ * * - ' (\s -r ,) u f 

the convergence for this series expansion is quite rapid and often values 



(6.9) 



of u and v of 2 will give good accuracy for short time lag and low frequenc- 
ies. 
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The Pade' approximation of e for u = v = 2 and if x = -sT is 



-Sr 


/ - §-(Sr) T~r, (5r^ 


(6.10) 


V ' — ' — v 


/ t $ CST)-> -3_ (ST) Z 


e ■ ST ^ 


- 6S7 + / 2 


(6.11) 




S*T 2 + 6 ST -t 12. 



the higher the powers of the numerator and denominator, the more terms of the 

series can be matched and the smaller the error in the time delay function. 

Equation 6,11 is the second order Pade 1 approximation. The fourth order Pade* 

- sT 

approximation (u = v = 4) for e is 

-sr (5J)* - 20 (sr) 3 + \$0($ T ) ~ 840 (sj) + /G80 

& ~ Z ; 5 (b . 12 ) 

{ST) + 20 (srf -r l8o( s ^)"4 S40 (sr) -t ! 680 

Higher order pade' approximations provide improved accuracy and fre- 
quency response. It is noted that the order of the characteristic equation 

- sT 

is raised by the order of the approximant used to simulate e . If the 
second order approximant is used, two additional roots are introduced into 
the characteristic equation. 

6.3 Simulation of transportation lag 
Equation 5.12 can be rewritten as 



_ ’20_ /_80 840 , I08o 

Jl-FT-V) = ' S7- (STf (srp _ CjrT (6vl3) 

~ ~ 2o J_ 8 o § 4 Q_ !_ 68 jO 

7 ST ' CSV* + (st)3 t (stJ* 

since the differential operator, D, and Laplace operator are interchangeable 
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when initial conditions are zero, equation 6.13 can be expressed as 



/ (* - T) 
-f(t) 



2_Q l_8o _ &4Q + t±£Q 
' D~f ' (.o ~) 2 (o r)* (or)* 



(6.14) 



/ + 20 + •8° 8±P + 1680 

' D 7 (or) 2 (OTp (Dr)*? 



rearrange equation 6.14 and we get 



i«-T)= -f(t) - + A^f(t-r)}-^ 




- f (t - rjJ -O. 






(6.15) 



according to equation 6.15, the analog computer circuit can be simulated 
as shown in Fig. 6.2, 

The time delay T in second is determined by the magnitude of the input 
resistors to the four integrators. Such a circuit works well only in a 
limited frequency range (direct current to approximately 50 cps). In- 
stability results if the variation of the input signal f(t) is too rapid. 

It is therefore not feasible to apply a step function to this circuit. In 
case f(t) is not equal to zero at the commencement of the computer run, it 
is necessary to apply appropriate initial conditions to the integrators. 

The transportation lag devices described thus far are limited to applica- 
tions in which the delay time T remains constant throughout the computer 
run. 

The method of simulating transportation lag by the Pade f approxima- 
tion yields what may be termed optimally flat phase characteristics at the 
origin but does not fit the phase characteristics as well over an extended 
frequency range. The major drawback to the use of the Pade 1 approximation 
to simulate transportation lag is their response to a step input or other 
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Pade approximac ion. 



input having high frequency components, since the amplitude ratio is 
always unity but the phase shift is approximated only for low C(f T pro- 
ducts. The phase shift characteristics of certain of the Pade' approximation 
are shown in Fig. 6.3. As a result of the inability of this approximation to 
linearly phase shift high frequencies, the output of an approximant rings 
badly for high-frequency inputs. This behavior is illustrated in Fig. 6.4. 

The actual location of the transportation lag transfer function in the 
forward path is immaterial and it may be placed at the end of the path in 
the computer where the high frequency components have been attenuated to the 
greatest possible extent. 

In the great majority of cases the fourth order Pade' approximant is 

more than adequate for the simulation of transportation lag. In cases 

where linear phase shift of higher ccrT products than those shown in Fig. 

6.5 is desirable, higher order Pade 1 approximants or more complex curve 

- sT 

fitting techniques may be used to simulate e . The phase shift character- 
istics of the sixth order and eighth order Pade 1 approximants are shown in 
Fig. 6.5. On the basis of this result it may be inferred that the order 
of the Pade' approximation may be increased to increase the magnitude of 
the cw T product that is linearly phase shifted without distorting the 
phase shift of low t products. Because the Pade 1 approximation may be 
expressed in a general form from which any order approximant may be develop- 
ed, it is felt that the use of curve fitting techniques requires a needless 
expenditure of effort without any tangible benefit. 

6.4 Generation of time delay from magnetic tape recorder. 

The transportation lag generated by using the magnetic tape recorder 
is very accurate, and it saves a lot of time and labor for simulation by 
using amplifiers. 
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WT Product (radians) 




Fig. 6-3 Phase Shift Characteristics of the Fade 

Approximat ion 
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(a) delay time ^ 1»0 second 




(b) delay time ^ 10 seconds 

Fig. 6-4 Curves obtained from Pade approximation 

l 
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Phase shift (radian) 




Fig. 6-5 Phase shif Oharacterist ics of the 
Pade Approximation 
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The magnetic tape recorder used has the following specifications: 
Ampex Model FR1107 magnetic tape recorder 
Tape speed: 3-3/4, 7-1/2, 15, 30, and 60 inches per second. 

Input impedance: 100,000 ohms unbalanced to ground . 

Input level: 1 volt rms recommended for normal recording level. 

Output impedance: 125 ohms 75 /yf, unbalanced to ground. 

Output level: 1 volt rms normal, across 10,000 ohms load impedance at 

normal recording level. 

The response curves obtained from the magnetic tape recorder are as 
shown in Fig. 6.6. Compare this with that obtained from analog computer 
as shown in Fig. 6.4, it is obvious that the magnetic tape recorder is 
most desirable for simulating the time delay, since it gives a more ac- 
curate result. 
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(b) delay time z. I*? seconds 
Fig. 6-6 Curves obtained from magnetic tape reoorder 

I 
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7 . System simulation . 



7.1. Analog computer simulation 

The block diagram of the system to be investigated is as shown in 

Fig. 7.1. 




Fig. 7.1. The block diagram of the relay servo system. 

The simulation of the transportation lag is shown in Fig. 6.2. The simula- 
tion of the ideal relay is as shown in Fig. 7.2. 




Fig. 7.2. The simulation of the ideal relay. 

From the block diagram of the system, the analog computer circuit 

diagram can be simulated as shown in Fig. 7.3. 

Because of its convenience and accuracy, the magnetic tape recorder is 

_ s t 

selected for the generation of transportation lag (e ) through the experi- 



ments. 
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Fig. 7-3 Operational Diagram of Analog Computer setup for the 
block diagram of the relay servo system of Fig. 7-1. 



7.2. Digital computer program 



The block diagram of the system is shown below. 



f(t) + _ + / - N 5(t) 






E (IT) 



Delay 




Relay 


b(t) ^ 


IDl 


EDELAY^ 


V 


B(IT) 




CTACH 



DA DELAY 



ADELAY 



DELAY 

ID2 



Fig. 7.4. The block diagram of relay servo system. 

The variable names used in the program, E(IT), EDELAY, etc., are shown 
below their mathematical counterparts. Note that the fixed point IT or the 
floating point T indicates the number of increments, starting with IT = 1 : 



zero time 

i 



r 



.time = DELAT time = 2 DELAT, etc. 



IT = 1 
T = 1. 
E (1 ) 



IT = 2 
T = 2. 
E (2) 




etc. etc. etc. 

The method of computation is simply to evaluate the analytic expression for 
the response at each increment point. The response is always supposed to 
have started at the time of the last relay operation (TF in increment measure) 
with initial condition equal to the final values of A(IT) = AF, and DA(IT)/ 
CTACH = DAF at the end of the previous response segment: 
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new response computation starts here 
based on relay operation at IT = TF 



j 



^if computation at this point shows change of sign of EDELAY from previous 
lue, the relay operates here, TF, AF, and DAF are recorded as the current 



va 



value of IT, A(IT) and DA (IT) . 
Analytic response. 

The plant transfer function is 

A(s) / 



B(V 



5(3 + <*) 



e(s)— S 2 /l(s) -t o(SA<s) 

and the differential equation describing the plant is 

b (t )= + <* 

with initial conditions a(0) and a(0) (AF and DAF in program). 
The input is a step V 

^ = 5 a /l(s) - Sa(o) - d(oJ + oCSACs) - o(a(o) 

A Vt S j q(Q) 

'• S 2 (S-h <x) 

This can be expressed into partial fractions 

A(s) = JL. + .JL _ X(s ■<■*) + vs(s-f«j+2s 2 

s 



5 S +<* 



S 2 (S ■+ ot ) 
2 



(7.1) 

(7.2) 

(7.3) 



(7.4) 



(7.5) 



equating the 


coefficient of s 


5^ CS-+ * ) 




S 2 : 


V -+■ z = 


a. (o) 


(7.6) 


5' : 


X ■+ Yoi 


d{o) -h o(CUo) 


(7*7) 


S° : 


x°< — 


V 


(7.8) 
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\/ v 

X = 2 . ^ - = V/ALPHA in program 

Y = 5 - jd(O) -f- o(.<Z(o) - )() = (DAF + ALPHA-X) /ALPHA 
^ = d(O) " Y = AF - Y 

from equation 7.5, the time solution is 

-at 

act)- xt + Y ze 

= X*TAU 4 Y 4 Z*EXPF ( -ALPHA *TAU) in program 

a(±)~ x -z<* 

= X - Z*ALPHA*EXPF(~ALPAH*TAU) in program 
t-^™^TAU is time from last relay operation 
Program accuracy. 

1) . The time delays can be introduced only as an integral number 
of increments (each of magnitude DELTAT) . 

2) . The only actual inaccuracy is the fact that the relay can only 
operate at integral values of IT*. The accuracy, for a given increment size, 
DELTAT could be fairly easily improved by introducing an interpolation pro- 
cedure between computation points. The same technique could be used to 
introduce non-integral number of the delays. 

The plant dynamics can be easily changed by re-working the partial 

fraction expression and the expression for a(t) and a(t). 

The program for the relay servo system is shown below. And one case 

of its output (t, = .8 second, k = 2,0) is attached to this program, 
d t 



*This error may be regarded as a variability in the forward path time delay, 
so that the actual delay may lie between the desired ID1 and ID1+1. 
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.JOB CHI ANG , MAX. TIME 3 MIN. 

PROGRAM TLRSOJ 

RESPONSE CALCULATIONS FOR RELAY SERVO SYSTEM WITH DELAYS ID1.ID2 AND 
ID3 IN THE FORWARD, FEEDBACK AND TACHOMETER PATHS. RESPECTIVELY. 
SOLUTION IS OBTAINED FOR A TOTAL TIME OF TMAX SECONDS IN INCREMENTS 
OF DELTAT SECONDS. IN THE PROGRAM BOTH T AND IT ARE THE NUMBER OF 
TIME INCREMENTS PLUS ONE. THAT IS, THE IN I TI ALr’ T IME IS T=1.. THE THREE 
DELAYS ARE MEASURED BY THE NUMBER OF TIME INCREMENTS. THE SYSTEM 
FORCING FUNCTION IS SPECIFIED BY PROGRAM STATEMENT NUMBER 1, AND TWO 
FOLLOWING CAROS. v 
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INPUT DATA. 

IN 

101 , 102 , 

V 

ALPHA 

CTACH 

TMAX 

OELTAT 

INCPRNT 

INCPRNT IS THE 
V IS THE RELAY 
ALPHA IS THE 
CTACH IS THE 



ON CARD 
ON CARDS 
ON CARD 
ON CARO 
ON CARD 
ON CARD 
ON CARO 
ON CARO 
NUMBER OF 
VOLTAGE 
PLANT TIME 
TACHOMETER 



1 

2 



RIGHT JUSTIFIED iM COLUMNS 1 AND 2 



THRU 4. RIGHT 



COLUMNS 

5. IN COLUMNS 1 THRU 10 WITH DECIMAL 

6. IN COLUMNS 1 THRU 10 WITH DECIMAL 

7. IN COLUMNS 1 THRU 10 WITH DECIMAL 

8. IN COLUMNS 1 THRU 10 WITH DECIMAL 

9. IN COLUMNS 1 THRU 10 WITH DECIMAL 

10. RIGHT JUSTIFIED IN COLUMNS 1 AND : 

INCREMENTS BETWEEN PRINT OUTS 



1 THRU 3 

POINT 

POINT 

POINT 

POINT 

POINT 



CONSTANT 

CONSTANT 



1 02, 1 03, V, ALPHA, CTACH, TMAX, CELT AT, INCPRNT 



THE MAXIMUM NUMBER OF INCREMENTS IS 1000. THE NUMBER OF INCREMENTS 
MAY NOT EXCEEO 900 IF A GRAPH IS REQUIRED. 

IN IS THE NUMBER OF CASES TO BE SOLVED. NOTE THAT THERE MUST BE ONE 
DATA OECK OF NINE CAROS FOR EACH CASE. THE FIRST CARD CONTAINING IN 
IS NOT REPEATED. PUT IN=1 IF ONLY ONE CASE IS TO BE SOLVED. 

EACH SET OF DATA CAROS MUST CONTAIN THE REQUIRED GRAPH FORMAT CAROS. 
REMOVE CALL GRAPH STATEMENT AND GRAPH FORMAT CARDS IF NO GRAPH IS 
REQUIRED. 

DIMENSION A{ 1001 > ,B( 1001 ) ,E( 1001 ) ,DA( 1001 ) ,DE( 1001 } ,ET ( 1001 ) , 

1 TIMES! 1001 ) ,FIN(1001 ) » 8S I 1001) 

READ 150, IN 
FORMAT (12) 

00 34 K = 1 , IN 
READ 100,101, 

PRINT 200 
PRINT 201 

PRINT 202, 101 ,102, 103, V, ALPHA, CTACH, TMAX, DELTAT 
FORMAT ( 1 3/13/ I3/F 10.1 /FI 0.1 /FI 0.1 /F 10.1 /F10. 1/12) 

FORMAT {//) 

FORMAT ( 3 1 H RELAY SERVO SYSTEM INPUT OATA //) 

FORMAT t 8H ID1 =13/ 8H 102 =13/ 8H 103 =I3/8H V =F7.3/ 

1 8H ALPHA =F7 .3/ 8H CTACH =F7.3/ 8H TMAX =F7.3,5H SEC./ 

2 ’ 8H DELTAT*F7. 3, 5H SEC.) 

PRINT 203 
PRINT 204 

FORMAT { 6X42HCALCULATED RESPONSE OF RELAY SERVO SYSTEM ) 

FORMAT (1H1) 

T = 1 . 

TF= 1 . 

A ( 1 ) =0. 

DAI 1 ) =0 . 

AF=0 • 

D AF=0 • 

EDELAY=0. 

INC=TMAX/DELTAT 
NUM= INC-2 
F= 1 . 

DF=0. 

FORMAT I6X30HF0RCING FUNCTION * A UNIT STEP) 

FINI IT)=F 
IFIT-1.) 7,2,7 
PRINT 205 
PRINT 200 
PRINT 206 
PRINT 207 

FORMAT ( 7X4HTI ME5X7HF0RC ING5X6HSYSTEM5X5HPLANT4X5HPLANT7X4HRATE 
1 7X4HRATE) 



150 



100 

200 

201 

202 



204 

203 



1 

205 



206 
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3 

4 

5 

6 



208 



207 FORMAT ( 1 6X8HFUNCT I0N4XSHERR0R6X5H I NPUT4X6H0UtPUT5X6H0UTPUT 
1 6X5HERROR/ ) 

E ( 1 ) =F 
ET ( 1 ) =F 
I F ( E ( 1 ) ) 3. 5. 3 
I F ( 101)5,4.5 
B ( 1 ) =SI GNF ( V, E( 1 ) ) 

GO TO 6 
B ( 1 ) *0. 

TIME*0. 

0E( 1 ) =DF 

PRINT 208, TIME,F,ET( 1 ), B( 1 ) , A ( 1 ),DA( 1 ),DE( 1 ) 

FORMAT (6XF6.3,3XF8.4,4XF8.4,3XF6.2,3XF8.4,3XF8.4,3XF8.4) 

DO 33 I T=2 » I NC 
T = I T 

TIME=(T-1 . )*OELTAT 
TIMES! I T ) =TI ME 
EDELAY1=EDELAY 
J = 1 

GO TO 1 

7 TAU= ( T-TF ) *DELTAT 
TAU IS TIME FROM LAST 
TF IS THE TIME OF THE 
X=B( IT-1 ) /ALPHA 
50 Y*(DAF+ALPHA*AF-X) /ALPHA 
Z=AF-Y 

AF AND DAF ARE THE PLANT OUTPUT AND RATE OUTPUT AT TIME OF LAST RELAY 
OPERATION. X, Y AND Z ARE THE COEFFICIENTS IN THE PARTIAL FRACTION 
EXPANSION. 

A( IT)=X*TAU+Y+Z*EXPF( —AL PH A* TAU ) 

DA (IT)* X-Z* ALPHA *EXPF (- ALPHA*TAU) 

I F ( 102)9.8,9 
ADELAY=A( I T ) 

GO TO 12 

IFUD2-IT+1 )10. 10, 11 
ADELAY=A( IT-ID2) 

GO TO 12 
ADEL AY=0. 

I F ( I D3 ) 14,13,14 
DADELAY=DA(IT) 

GO TO 17 

I F ( ID3- IT+1) 15*15,16 
DADELAY=DA( IT-ID3J 
GO TO 17 
D ADELAY*0. 

E( IT )=F-ADELAY-DADELAY*CTACH 



RELAY OPERATION (OR FROM TIME ZERO) 

LAST RELAY OPERATION (INCREMENT MEASUREMENT) 



8 

9 

10 

1 1 
12 

13 

14 

15 

16 
17 



51 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 
210 

29 

30 

331 

332 

333 



GO TO (51,31 ),J 
I F ( ID1 ) 19, 18,19 



EDELAY=E( IT) 

GO TO 22 

I F ( ID1-IT+1 ) 20.20, 21 
EDELAY=E( IT-ID I ) 

GO TO 22 
EDEL AY=0. 

IF ( E DEL AY 124,23,24 
B ( I T ) =0 . 

GO TO 25 

B(IT)=SIGNF(V,EDELAY) 

IF ( XMODF ( ( IT-1 ),40*INCPRNT> ) 27, 26,27 
PRINT 203 

IF (XMODF ( ( IT-1 ) , 10*INCPRNT) ) 29,28,2 9 
PRINT 210 
FORMAT (/) 

I F(8( I T ) — B ( IT-1 ) ) 30, 31 *3 0 
IFlIDl-IT+l) 332,331 , 332 
S=0. 

GO TO 333 

S=ABSF(EDELAY)/(ABSF(EDELAY)+ABSF(EDELAY1) ) 
TF=T-S 

AF*A(IT)-S*(A( IT )— A( I T— 1 ) ) 



I 
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DAF = DA( IT)-S*(DAUT)-DA{ IT-1)) 

C THE LAST SEVEN STATEMENTS SET THE INITIAL VALUES FOR THE NEXT SEGMENT 
C OF THE CALCULATION IF RELAY OPERATION HAS BEEN DETECTED BY STATEMENT 
C NUMBER 29. THE TIME OF SWITCHING HAS BEEN FOUND ‘BY LINEAR 
C INTERPOLATION AS HAVE AF AND DAF. F SHOULD BE CONTINUOUS EXCEPT AT 
C TIME ZERO. '* 

T IME=( T-S-l • ) *DELTAT 
PRINT 209. B( IT) .TIME 

209 FORMAT (22H THE RELAY SWITCHED TO F6.1.10H AT TIME = F7.3.5H SEC.) 
C WE HAVE NOW TO CORRECT A ( IT) . 8( IT ) AND EC IT) FOR SWITCHING AT TF. 

T IME= (T- 1 • )*OELTAT 
J=2 

TAU= (T-TF)*DELTAT 
X = B ( IT ) /ALPHA 
GO TO 50 

BSC IT)=B( IT)«0. 125/ABSFC V)+0. 125 
ETC IT)=F-A( IT) 

DEC IT)=DF-DA(IT) 

IFCXMODFC ( IT-1 ) , INCPRNT) ) 33.320. 33 
PRINT 208, T I ME, F, ETC IT), BC IT ) , AC I T ) ,0A( I T ) , DEC I T ) 

CONTINUE 

( NUM.ET ,DE,8) 

( NUM, A , DA , 8 ) 

(NUM.TIMES.A, 8) 

(NUM, TIMES, DA, 8) 

(NUM, TIMES, BS, 8) 

(NUM, TIMES, ET, 8) 
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320 

33 



34 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

PRINT 

END 

END 



GRAPH 

GRAPH 
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7.3 Comparison of the results obtained from the digital computer and from 
the analog computer. 

The result obtained from the digital computer for the system with time 

delay .4 second has been checked with analog computer. They are quite in 

agreement with each other as shown in Fig. 7.5 and 7.6. 

comparison of the results 

by digital by analog 

computer computer 



overshoot and undershoot 3.5 2.6 

rise time 5.9 5.7 

Because the digital computer can give an accurate and exact value for 
any desired instant, therefore, the digital computer is selected for con- 
ducting most of the experiment work, and the analog computer is used to 
check the results obtained by digital computer. They are quite in agreement. 
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Ob) results obtained from analog computer 
Fig. 7-5 Output, velocity of the system, t^ - . 
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4 second 
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8. Effect of transportation lag on relay servo system . 



8.1. General 

From the proceeding sections* we know that the time delay has an adverse 
effect on control systems* and which is unavoidable in nearly all the control 
systems. Because the time delay may basically be represented by an infinite 
series of terms* this creates an infinite number of roots. 

Later on in this section, it will be noted that the time lag results in 
successive overshooting and undershooting of the correct final value* re~ 
presenting serious oscillations about the final position. In addition to 
this it will also change the root locus and phase trajectories. 

8.2 The frequency response. 



From Fig. 4.6(c) we have the transfer function with removed 



Gc y 



f / -t- /C t 5 .) g 
S (S ■+ ■ I ) 



-ST 



( 8 . 1 ) 



where K_ and T are variables. In case of time delay equal to .5 second* it 
is found that the tachometer gain, = 1.5, will compensate the system to 
an acceptable optimum operation, this will be seen’ in section 9, now equation 
8.1 can be written as 



G(s)= 

5(S-+‘) 



( 8 . 2 ) 



The Bode diagram for the above transfer function is shown in Fig. 8.1. So 
far as the linear system is concerned* the gain margin for .5 second time 
delay system is approximately 1.6 db, but for the undelayed system* it is 
stable for all gains and never unstable. 

The time delay has introduced additional branches into the diagram. 
These branches appear due to the periodic nature of the time delay* and the 
increasing phase shift introduced at higher frequencies. It is evident that 
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Fig* 8-1 Bode diagram of Eq c 



the gain margin for each of the branches increases* This is analogous to 
the statement that the roots of the primary branch would go unstable first 
if the gain were raised* 

The gain vs* logos curve of the delayed system is identical to that 
of the undelayed system, since the time delay introduced no additional gain 
to the system* The phase angle vs* log^u* curve of the delayed system is 
drastically changed since the equation for the phase angle is 

-( 90 ° -h OJ -f cot) ~t~ tan' (/■ 5co) (8.3) 



The CO T term in the above equation has a predominate effect to the phase 
angle, it introduces additional phase shift and extra branches* 

The limit cycle occurs at OJ = 1.7, k = 1*20* 

8,3 The root locus. 

The root locus of any system will be greatly changed when a transporta- 
tion lag was introduced. As a matter of fact, for an undelayed system the 
transfer function in Equation(8. 2) becomes 

_ A 5 ( s ' & y) 



C? ( $) 



s CS-hi) 



(8*4) 



no limit cycle exists, because the root locus never intersects the imaginary 
axis as shown in Fig* 8.2, which shows a stable system for all gains* 




But when time delay is introduced this situation will be changed* The root 
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locus for Eq, 8,2 (,5 second delayed system) is shown in Fig. 8,3, The 
construction details for this root locus were described in section 4,4,2. 

Fig, 8,3 illustrates the effect of time delay on the root locus. It 
is apparent that this system will go unstable with only a small increase 
in gain. The primary branch stability limit before the root locus corsses 
the imaginary axis are k = 1,2, , and OJ - 1,76 rad/sec, , where the limit 
cycle occurs, this almost checks with Bode diagram where k = 1,20, <xr = 1,7 
rad/sec. 

Because time delay introduces periodically phase shift of 360 degrees, 
so there presents additional (upper) branches when at high frequencies. It 
is noted that for the same gain the additional roots on the upper branches 
are attenuated more than the roots on the primary branch. This substantiates 
the Bode diagram findings, which shows that the roots on the primary branch 
phase shift would cause instability, 

8.4 Transient response. 

Because of time delay introduced in the system, there is no graphical 

method that can be easily used to produce the transient response of any 

system. The transient response of the delayed relay servo system (second 

order) for various delay time and various k (tachometei gain) are shown 

t 

in Fig. 8.4. In Fig. 8.4(a) there are six curves under the same taebo- 
meter gain (k fc - .02); curve (1) is an undelayed system, the rest of the 
five curves are with ,1, .25, ,5, ,8 and 1,0 second delay time respectively. 
For 1 second time delay system it presents an overshoot and undershoot of 
approximately 60 percent. 

Table 8,1 is the data obtained for M , t , t under different values 

P t P r 

of and t^, Plotting these data into curves it is found that they appear 

as a straight line. In Fig. 8.5 are shown M and t versus delay time, in 

P t P 
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Fig. 8-3 Root locus of relay servo system with time 

delay of .5 second. 
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-4 Cl, Response of relay servo system t,zO, •!» .25, .5, .8, 1.0 second respectively. 
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Fig. 8 -4b , Response of relay servo system with t d= 0, .1, .25, .5, .8, 1.0 second 

respectively, k. -.8 




Table 8,1. Data for first overshoot t and M under various time delays 

p p t 

for a step input. 



K t 


T d 


M 

P t 


t 

P 


t 

r 


amplitude of 
limit cycle 


.02 


.0 


1.1285 


3.50 


2.828 


no limit cycle 


If 


.1 


1.1870 


3.70 


2.95 


.08 


II 


.25 


1.2510 


4.10 


3.095 


.1625 


II 


.5 


1.3750 


4.60 


3.350 


.2575 


II 


.8 


1.5300 


5.10 


3.650 


.4910 


II 


1.0 


1.6500 
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t (time for first overshoot) 
P 



- - - - M (magnitude for first 



overshoo t.) 
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Fig, 8,6 is shown the rise time versus delay time, there it shows that the 
effect of time delay on the transient response is to increase the over- 
shoot and undershoot and cause continuously oscillations and make the rise 
time longer, it turns out that the response of the system can never be 
settled. 

It is noted that M , t , t increase as the delay time increases, 

P t P r 

But when we increase the value of k^_ they decrease. And to a certain value 
of k some of them become acceptable optimum operation this will be seen in 
section 9.2. 

In Fig, 8.7 are shown the plant output, rate output and system error. 
Fig. 8.7(a) is the response of the undelayed system, the plant output set- 
tled at 7 seconds; Fig. 8.7(b) shows the plant output, rate output and 
system error for a system with delay time of ,5 second and k^ = ,02, the 
output continuously oscillates with an overshoot and undershoot of approxi- 
mately 30 percent; Fig, 8.7(c) shows the plant output rate output and system 
error for a system with time delay of one second and k^ - .02, The plant 
output continuously oscillates with an overshoot and undershoot of approxi- 
mately 60 percent, the frequency of oscillation is slower than that of ,5 
second time delayed system; Fig, 8.7(d) shows the plant output, rate out- 
put and system error for a system with delay time of 2,5 seconds and k^_ = 
.02, the plant output continuously oscillates with an overshoot and under- 
shoot of approximately 135 percent, the frequency of oscillations is much 
slower than that of .5 second time delayed system, 

8.5 Time delay in the phase plane diagram. 

The effect of time delay in application or removal of the applied 
torque after the control has signaled that the operation should be perform- 
ed. It will be of interest to study the effect in terms of the phase plane 
diagram. 
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Fife. 8-6_ J Time rise vs. delay time 
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plant output, velocity /S nd system error of 'an undeleyed relay servo system 
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Fig. 8-7 b • The plant output, velocity end error of relay servo system vdth a time delay of 
second and k - .02. 
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plant output, velocity and system error of a relay servo system with a tifl9 delay 
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Because of the time delay the trajectory in the phase space will be 
prolonged beyond the point c° = R* by an amount such that the difference in 
coordinates is just equal to the time delay* that in Fig* 8*8 the torque 

should reverse at a s but it does not reverse until B is reached* and (C^ * C^) 
(v® - v')/2 is equal to t® * Similarly*, the torque is reversed at point D 

D A Q 

instead of at c 5 and at F instead of at E* It is apparent that a curve could 
be drawn through points such as B and F so that the switching of torque 
direction would take place when the trajectory crossed this line 8 and the 
line from B to D would play the same role as the line through R° for the 
control without time delay* 




Fig* 8*8. Normalized phase plane diagram for relay servo 
system with a time delay of *1 seconds for a 
step unit* 

The effect of time delay is thus to skew the switching line from the 

verticle in the phase space s and quite approximately the angle of skew is 

-tan *t , and is in clockwise direction* 
a 

It is seen in Fig. 8*8 that the trajectory appears to repeat itself^ 
since the last two oscillations are nearly identical* The control oscillates 
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with a constant amplitude on the trajectory. Such a repet ive diagram re- 
presenting the steady state oscillation is called limit cycle. From the 



experimental work it shows that the relay controls with time delay and nr. 
dead zone exhibits the defect of continuous oscillate loins . 

In Fig* 8.9 is shown a family of trajectories so arranged that the 
initial and final points of each are symmetrically located with respect to 
the origin. Thus each trajectory represents just such a limit cycle* and 
since the time delay for each is fixed* the locus of the end points of these 
limit cycles measures the time delay required to obtain the amplitude of 
oscillation. As would be expected* the larger the time delay the larger is 
the magnitude of the oscillation* and with that* the lower is the frequency 
of oscillation as indicated in Fig, 8„5„ In Fig. 8.10 is shown a plot of 
time delay and oscillation amplitude 0 It turns out to be a straight line* 
therefore, ve can obtain an equation to express the relation of oscillation 
amplitude with time delay as 



It Is interesting to see the relay type control that can continuosuly 
cycle when a time delay is present. The relay type of control always 
oscillates about the equilibrium position. For any time delay the system 
would always oscillate* and with a magnitude of oscillation depending on 
the time delay. The relay type control can never be unstable in the sense 
that the amplitude of oscillation will continue to increase without limit. 
It is obvious that this is the advantage of the relay servo system. 

In Fig, 8.11(a) the phase plane of the relay servo system with a time 
delay of zero and .5 seconds is plotted on the same sheet. In Fig, 8, 11(b) 



Aosc * & 



(8,4) 



the above equation is for a second order system specified by 




J 
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Fig.'S - 9 . Phase plane trajectory showing loci of limit cycle for oscillation of constant 
amplitude with various time delay 
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Fig.' 8-lla • Phase plane diagram for e relay servo ssystem showing the effect of time delay 





Fig. E-lib. Phase plane diagram for a relay servo system showing the effect of time 
delay 






is the phase plane of the relay servo system a time delay of one and 2*5 
seconds plotted on the same sheet. From these phase portraits we find a 
most important aspect of the phase portrait is the switching line which is 
not vertical. The effect of the dead time is to rotate the switching line 
into the first and third quadrants* the result is growing oscillation* 

The process of moving the switching line from vertical axis to a line 
in the first and third quadrants is seen to affect the nature of the response 
very unfavorably. On the other hand* if the switching line can be made to 
lie in the second and fourth quadrants by means of phase advance* or lead., 
on the error signal* the relay servo system is likely to be a very well be- 
haved one* in that it has small overshoot and will reach and remain within 
a given value of error in a short time. This will be seen in section 9* 

Dead time will result in a terminal limit cycle which may or may not be 
acceptably small* the amplitude of oscillation depends upon delay time as 
shown in Fig. 8.9 and .10* 

8.6. The effect of time delay in steady state oscillation frequency. 

In Fig. 8.12 the curve of system error against time is shown* Initial 
condition of the error is due to a suddenly applied step input signal* this 
gives an initial error and velocity (output rate) to the form of response 
waveform as shown in Fig. 8.12. After the transient has died out the out- 
put shaft follows the input by performing oscillation about a final posi- 
tion as shown in Fig. 8.5. 

In Fig. 8.12 curves 1* 2* 3, 4* and 5 are for a delay time of .1* -25* 

.5, .8 and 1.0 second respectively. It indicates that the smaller time 
delay causes high frequency oscillation with smaller amplitude* in other 
words* the frequency of the oscillation is higher* and the larger value 
of time delay results in larger oscillation amplitude with lower frequency 
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as shown in Fig, 8.12, 

In Fig. 8.13 curves 1, 2, 3, 4, 5, and 6 are for a delay time of 0> . 1, 
.25* ,5, .8, and 1.0 second respectively. The curves show the relation of 
the rate output to time delay. 

From Fig. 8.5 or Fig. 8,12 S the data for oscillation frequency against 
time delay can be obtained as shown in Table 8.2. 

Table 8.2, Fata for steady state oscillation frequency vs. time delay. 



delay time 

.1 

.25 

.5 

.8 

1.0 



frequency 

,455 

,294 

.200 

.156 

.135 



The data plotted in log log scale sheet there appears to be a straight 
line as shown in Fig. 8,14; it is seen that the frequency decreases as 
the delay time increases. Because the curve of oscillation frequency 
against time delay is a straight line, therefore we can obtain an equation 
to express their relation as 



*/3S 

f~ - s7 (8.5) 

L d 

From Fig. 8.12 we take the steady state error for each time delay, 
tabulate the data as shown in Table 8.3. 




Table 8.3. Steady state error vs. time delay. 



Delay time Error 

0 0 
.1 .09 

. 25 . 18 

.5 .31 

.8 .484 

1.0 .588 



Fig. 8.15 is the plot of error vs. delay time. It comes out approximately 
a straight line, it indicates that the steady state error is proportional 
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l2JSrror vs. time curves showing the effect of time delay on the system error 
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Figj 8 ™ 13 c Velocity vs. time curves showing the effect of time delay on velocity. 





Fig 8-14*. Steady state oscillation 
frequency vs. time delay. 




Fig. '8-15? Systemerror against time delay 



1.0 



Co the t i iK«ii delay. 

8.7. The effect of the magnitude of the relay voltage. 

From the experimental work, we find chat the magnitude of the relay 
voltage also has some effects on the response curve and phase trajectory of 
the system. VJhen the relay voltage is smaller the response has a longer 
rise time and a smaller overshoot and undershoot, and there is presented 
some velocity saturation in the phase trajectory* the amplitude of the 
limit cycle is smaller. 

In Fig. 8,16 shows two curves with the same delay time of *5 second 
and same value of k, - ,02* but with different relay voltage, curve (1) is 
with V volt and curve (2) is with V ~«5 volt. From these curves it is 
indicated that the smaller relay voltage causes smaller oscillation amplitude 
and lower frequency, but it presents a longer time rise which is its dis- 
advantage; so far as the time rise is concerned the larger value of relay 
voltage is desirable, therefore, in all the experiments we set the relay 
voltage equal to .5 volts. The rise time for .2 volts of relay output is 
6.5 seconds, and the rise time for .5 volts is 3.35 seconds. The overshoot 
and undershoot for .2 volt system is 12,5 percent, and that for .5 volt is 
30 percent*. 

In Fig. 8.17 the phase plane with the same delay time and same value of 
is shown, but with different relay voltage, it presents a velocity satura- 
tion for the system with .2 volts of relay output, and the limit cycle is 
smaller. 
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16. Curves showing the effect o£ the relay voltage to the response. t. = .5,k^=.02 
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Fig. 8-I?.Fhase plane diagram showing the effect of relay voltage, t.r.5 second, k t= . 02. 




9 . Elimination of the steady state oscillation 



In many applications steady state oscillations of the output shaft are 
inadmissible; besides causing continual operation of relay contacts and 
wear of gears and bearings the power consumed in oscillations may be objec- 
tionable. The simplest method of eliminating the oscillation is to intro- 
duce a dead zone in the relays so that the error always returns to the 
dead zone after a disturbance. But the relay used has been thus far assumed 
an ideal relay with no dead zone, so we shall not discuss this possibility, 

9.1 Velocity feedback stabilization 

A very effective method of reducing the settling time and limit cycle 
oscillation of an relay-servo is to make the relay switching conditional 
not only on misalighment but also on its rate change. This can be achieved 
by feedback within the position control loop if the feedback signals can 
be added algebraically to the misalignment signal at the error sensing 
device . 

The effect of negative velocity feedback is illustrated on the phase 
plane diagram in Fig. 9.1; and the corresponding time response is shown 
in Fig. 9.2. 

It is obvious that the larger the value of k the smaller the ampli- 
tude of the limit cycle, and the more the switching line shifted to counter- 
clockwise direction; and the amplitude of oscillation is reduced; the over- 
shoot and undershoot for k^_ = .02 is about 15 percent, for k^ = .8 is 
about 2 percent as shown in Fig. 9.2, 

9.2 Optimising the transient response. 

Because of the introduction of the time delay in the system, it is 
difficult to derive a method to determine the value of k^ which will give 



101 




102 



Fig. 9 = 1 0 Phase plane of relay servo system with a time delay, of .25 second compensated 
by increasing the gain of tachometer feedback gain to .8 
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Fig. 3 - 2 . Curve showing the response of the relay servo system be compensated to an acceptable 
optimum operation by in creasing k t to a value of .8, time delay .25 second. 



an optimum operation, which can only be determined through experimental ion* 
As we know the system damping will be improved by the tachometer feed- 
back, the net effect of the derivative signal is to cause the relay to 
operate early, so by using k^ we can control the switching characteristics 
of the relay. For many applications of relay servomechanism, the tacho- 
meter feedback compensation is adequate* 

Through experimentation it was found that the system with a delay time 
below ,5 second can be compensated for 5 percent overshoot and undershoot 
which is considered to be acceptable. 

The response and phase plane of the compensated system for various time 
delays are shown in Fig* 9.3(a) through Fig* 9. 3(e)* They have been compen- 
sated by velocity feedback. The values of k (which give acceptable optimum 

t 

operation) are as follows: 

Table 9*1 Optimum value of k for eliminating the effect 

of time delay* 

t k o/o overshoot 

and undershoot 



0 


.3 


/ 


.1 


.4 


.1 


.25 


.8 


1.6 


.5 


1.5 


5.2 



For systems with .8 and 1.0 second time delay , it is found that they can 

not be compensated for 5 percent overshoot and undershoot, even when we 

increase the tachometer gain, , to a considerable value, the amplitude 

of steady state oscillation is still rather large as shown in Fig. 9*4 and 

Fig. 9*5; the only effect of increasing k at this time is to shift the 

t 

whole response curve downward, as shown in Fig* 9*5 the response of the 

system with one second time delay at k - 4*5 is rather lew but still with 

t 
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considerable oscillation amplitude which is almost the same as that when 
= 1.5, which is beyond the acceptable value. 

Fig. 9.6 illustrates the time delay vs, curve at the range from 0 to 
.5 second time delay. The curve appears to be a straight line, but when 
the time delay above ,5 second occurs this situation does not hold, since 
the system with a delay time above .5 second cannot be compensated. When 
the time delay is less than .5 seconds, the system can be compensated; and 
because t^ vs. k^ curve is a straight line, therefore, it can be expressed 
by an equation as given in equation 9.1. 

K t = -+ ^-8 t d (9.1) 

Table 9.2 and 9.3 are the results for comparing the uncompensated and 
compensated system. The data in Table 9.2 is for the undelayed system 
for which k^ = . 3 will compensate the system to give optimum operation. 

The data in Table 9.3 is for a .5 second time delayed system for which k^ 
=1.5 will compensate the system to give acceptable optimum operation. 



105 




plant output, veloci ty and error vs. time 




phase plane (E vs E) 



Fig. 9 3a. Response and phase plane, 0, k^. - .3 
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plant output , veloci ty and error vs. time 




Phase plane 

Fig/'9-3b."J_ Response and phase plane, t^ - .1 sec. 
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phase . plane 

Fig. '9* 3c. ' Response and phase plane, - .25 sec. k^. - .8 
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plant output, velocity and error vs. time 




Fig. REsponse and phase plane, t^- .5 sec. - 1.5 
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1C. The response of relay servo system to ramp function input. 



In the previous sections we have been concerned primarily with the 
response to step function input. Experimental results obtained from the 
relay servo system with a ramp function input are shown in several graphs. 

When a unit ramp input is applied to the relay servo system, the out- 
put and the error of the system are shown as a straight line, and the velocity 
for most part is a horizontal line, as shown in Fig. 10.1. The Cdot vs. C 

and Edot vs. E phase planes have the same shape, merely the Cdot vs. C plane 

appears above X- axis. Fig. 10.2 shows the Edot vs. E plane. This is with 
a relay voltage equal to .5 volt. When the relay voltage is increased to 
1.0 volt the same situation remains. For relay voltage at .5 and 1.0 volt 

the only effect of the time delay on the relay servo system is to cause some 

delay in the response, no oscillation has been observed. 

But when the relay voltage is increased to 1.5 volts, this situation is 

* 

changed. For an undelayed system the output still appears as a straight line 
but the system error is no more a straight line, and the velocity curve is 
with some oscillations as shown in Fig. 10.3. Cdot vs. C and Edot vs. E 
phase planes have a different shape as shown in Fig. 10.4 and Fig. 10.5. 

For delayed system the output no more remains a straight line but with some 
oscillations, and the system error and velocity are much different than that 
of the system with relay voltage of .5 or 1.0 volt as shown in Fig. 10.6. 

The Cdot vs. C and Edot vs. E phase planes are quite different as shown in 
Fig. 10.7 and Fig. 10.8 respectively. Oscillations and limit cycle has been 
observed when the relay voltage is 1,5 volts. 

The effect of time delay on the response of the system can also be com- 
pensated by velocity feedback. For the .5 second time delayed system, when 
k = 2.0, the response can be compensated to an acceptable optimum operation 
as shown in Figures 10.9,10, II. 
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Fig. 10-1. C, CD0T,E vs. time 
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Fig. 1 . Phase trajectory of relay servo system, relay voltage 

.5 volt t- 
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Fig 15 - 7 . phase plane (rate output vs. output) of relay servo system with unit ramp input, 
time delay .5 second, k .02, relay voltage/. 5 volt. 




Fig. !£■ Phase trajectory of ralay servo system to an unit 
ramp input, time deiay — .5 second, relay voltage—/. 5 volt. 
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11. Conclusions, 



The aim of this thesis was to demonstrate the effects of transportation 
lag on the relay servo system. 

For conducting the experiment by means of an analog computer* first one 
could simulate the time delay by Pade*s approximation or by magnetic tape 
recorder. It was seen that the generation of time delay from magnetic tape 
recorder presents a more accurate and better result than that from the Pad e B s 
approximation. Therefore, the magnetic tape recorder was selected to 
generate the time delay throughout the experiment. 

In this investigation most of the experiments were conducted by using 
a digital computer which gave more accurate results than the analog com- 
puter does; especially since it can give an exact value at any desired in- 
stant. The analog computer has also been used to check the results ob- 
tained from the digital computer. It is noticed that they are quite in 
agreement . 

From the analysis and experiments the following conclusions has been 
obtained. 

The root locus of the system was greatly changed due to the intro- 
duction of the transportation lag; and because time delay introduced 
periodicaly a phase shift of 360 degrees, it results in an additional (upper) 
branches when at high frequency. 

The effect of time delay on the transient response is to increase the 
overshoot and undershoot and causes oscillations and also leads to a longer 
rise time, the response of the system can never be settled. 

The effect of time delay on the phase plane is to rotate the switching 
line into the first and third quadrants, the results is growing oscillations 
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In other words the transportation lag will result in a terminal limit 
cycle which may or may not be acceptably small. 

The plot of time delay versus oscillation amplitude came out to be a 
straight line; the delay time multiplied by a constant of .6 gives the 
magnitude of the oscillation amplitude. 

The time delay versus steady state oscillation frequency when plotted 
on log log coordinates came out to be a straight line in such a relation 
that a constant of ,135 divided by gives the oscillation frequency, 

A 

The steady state error and the velocity are proportional to the time 
delay. 

These effects caused by the transportation lag can be compensated to 
some extent. The method used in this investigation to eliminate steady 
state oscillation and to optimise the transient response by employing tacho- 
meter feedback. It was seen that a system with a delay time below ,5 second 
can be compensated to an acceptable small oscillation amplitude (5 percent). 
But for time delay above .8 second it cannot be compensated to such a small 
amount; the only effect of increasing tachometer gain is to shift the whole 
response curve downward with the same amount of steady state over- and under- 
shoot . 

The system with a time delay below .5 second can be compensated to reach 
an acceptable optimum operation. For delay time below 0 5 second its relation 
with tachometer gain is approximately a straight line, therefore, it can be 
expressed by an equation. 

When a time delay is present the relay type control can continuously 
oscillate, and with a magnitude of oscillation depending on the time delay. 
The relay type control can never be unstable in the sense that the amplitude 
of oscillation will continue to increase without limit. This is one advan- 
tage of the relay servo system. 
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